Bose Gases Near Resonance: 
renormalized interactions in a condensate 
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Bose gases at large scattering lengths or beyond the usual dilute limit for long have been one of 
the most challenging problems in many-body physics. In this article, we investigate the fundamen- 
tal properties of a near-resonance Bose gas and illustrate that three-dimensional Bose gases become 
nearly fermionized near resonance when the chemical potential as a function of scattering lengths 
reaches a maximum and the atomic condensates lose meta-stability. The maximum and accompa- 
nied instability are shown to be a precursor of the sign change of g%, the renormalized two-body 
interaction between condensed atoms. That is <?2 changes from effectively repulsive to attractive 
when approaching resonance from the molecular side, even though the scattering length is still pos- 
itive. This occurs when dimers, under the influence of condensates, emerge at zero energy in the 
atomic gases at a finite positive scattering length. We investigate the properties of near-resonance 
Bose gases via applying a self-consistent renormalization group equation which is further subject 
to a thermodynamic boundary condition. We also comment on the relation between the approach 
here and the diagrammatic calculation in an early article [Phys. Rev. A 85, 023620(2012)]. 



INTRODUCTION 

In the fascinating many-body world, the interplay be- 
tween few- and many-body physics is known to play a 
very important role in a variety of systems and is in- 
strumental to our understanding of diversified emergent 
phenomena. The best known example perhaps is the 
superconductivity phenomenon where the Cooper bound 
states due to the Pauli blocking effect lead to pairing, cor- 
relations and therefore superconductivity in solids [l], Q ■ 
Quantum degenerate gases of ultra cold atoms near Fes- 
hbach resonance that have been explored in many labo- 
ratories during the past decade [3[ are another fantastic 
platform to scrutinize the intriguing interplay between 
few- and many-body physics. Although far away from 
resonances, the physical properties of quantum gases can 
be well described by the existing dilute gas theories, very 
little, except in a few limiting cases, has been known 
about the fundamental properties of quantum gases near 
resonance beyond the usual dilute limit. 

The complexity of this topic to a large extent appears 
to be two-fold: 

1) the role of few-body physics (two-body, three-body 
states etc) in many-body systems. In other words, to 
what extent does the underlying few-body physics in- 
fluence the many-body correlations and which channel 
dictates the many-body properties near resonance? 

2) the effect of many-body background on the few-body 
structures. How are the few-body structures or multiple 
scatterings affected by the presence of many other iden- 
tical particles? 

These two issues can be treated successfully and sep- 
arately in the usual dilute limit. For instance, in the 
leading order, one can neglect the effect of finite-density 
background atoms on multiple scatterings or underlying 
few-body structures; the energy density can therefore be 



calculated perturbatively by assuming the few-body scat- 
terings are given by their properties in the vacuum and 
applying the low density expansion. The effects of few- 
body structures on many-body physics can be explored 
perturbatively. One can also further study the leading 
effect of quantum gases on dimers and trimers because 
the many-body states in the dilute limit are well-known. 
However, near resonance, these two issues are generi- 
cally entangled and ideally have to be addressed self- 
consistcntly, which usually becomes very challenging. 

In this article, we make an attempt to understand the 
fundamental properties of Bose gases near Feshbach reso- 
nance via examining the intriguing interplay between the 
few- and many-body physics in Bose gases at large pos- 
itive scattering lengths. For this purpose, we introduce 
a simple self-consistent renormalization-group-equation 
approach to address both sides of the coin. Many-body 
properties of a quantum gas are shown to influence the 
renormalization flow of few-body running coupling con- 
stants resulting in the change of the sign of the effective 
two-body interaction constants. That in return com- 
pletely dictates the many-body physics near resonance 
and leads to peculiar features in the chemical potentials. 
And we limit ourselves to the resonances with a very 
short effective range. 

A key concept we are going to focus on is the ef- 
fective interaction between condensed atoms near reso- 
nance. The common believe is although the underly- 
ing short range resonance interaction has to be attrac- 
tive, at low energy scales the two-body interaction is ef- 
fectively repuslive when the scattering length a is pos- 
itive. The argument runs as the following. The phase 
shift due to a short range attractive potential is given as 
5 = — arctanfca which yields 5 = —ka when k is small. 
The phase shift approaches — 7r/2 when ka is much bigger 
than unity but much smaller compared to 1/R*, R* (<SC a) 
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is the range of the attractive interaction. The low energy 
phase shift 6 — — ka turns out also to be the phase shift 
of any repulsive interaction with the same positive scat- 
tering length. And so the atoms interact effectively re- 
pulsively if one is interested in the scatterings at small k. 
However, the phase shifts of a short range attractive po- 
tential and a repuslive interation can become significantly 
different whenever ka is bigger than unity; for instance, 
for a hardcore potential with radius R — a, S — —ka for 
all momenta, differing from the value of —n/2 for aft- 
tractive potentials when ka ^> 1. Near resonance, a ap- 
proaches infinity and the issue of the effective interaction 
between condensed atoms becomes very subtle. Indeed, 
near resonance we show that condensed atoms no longer 
interact with an effective repulsive interaction even when 
the scattering lengths are positive, a somewhat surprising 
conculsion to many. 

The theory of dilute Bose gases has a long history. The 
theory for weakly interacting bosons was established by 
Bogoliubov more than half century ago 0, HJ . A properly 
regularized theory of dilute gases of bosons with contact 
interactions was first put forward by Lee, Huang and 
Yang [f| and later by Beliaev Q who developed a field- 
theoretical approach. Higher order corrections were fur- 
ther examined in later years Since these results 
were obtained by applying an expansion in terms of the 



small parameter 



(here n is the density and a is 



the scattering length), it is not surprising that formally 
speaking each of the terms appearing in the dilute gas 
theory diverges when the scattering lengths are extrapo- 
lated to infinity. This aspect of the dilute gas theory, to a 
large extent, is the main reason why a qualitative under- 
standing of Bose gases near resonance had been missing 
for long. In fact, the issue of quantum gases beyond the 
usual dilute limit has been one of the most outstanding 
unresolved problems in the quantum many-body physics. 

A few preliminary experimental attempts have been 
recently made to explore the properties of Bose gases be- 
yond the dilute limit thanks to the availability of Fesh- 
bach resonances in many cold-atom laboratories Ill4l4|. 
It has been suggested that when approaching resonance 
from the molecular side, Bose atoms appear to be ther- 
malized within a reasonably short time, well before the 
recombination processes set in, so to form a quasi-static 
condensate. Furthermore, the life time due to the recom- 
bination processes is much longer than the many-body 
time scale set by the degeneracy temperature which mo- 
tivates us to explore the thermodynamic properties of 
the quasi equilibrium gases. Here we intend to illustrate 
the intimate entanglement between few-body structures 
and many-body physics in a single approach. These il- 
lustrations suggest that most of the surprising features of 
near-resonance Bose gases can be understood in terms of 
the spectrum flow of few-body structures such as dimers, 
trimers which can be strongly influenced by the inter- 
actions with the condensates. The device we rely on to 



draw a qualitative picture is the renormalized interaction 
constant. 

The two-fold complexities mentioned above are actu- 
ally two sides of one coin and thus can not be com- 
pletely separated, although many theoretical attempts 
towards better understanding of near-resonance unitary 
gases have been focused on either of these two aspects. 
For instance, the role of three-body physics in Bose gases 
in the dilute limit has been explored and the contribu- 
tion of Efimov states 15-1|| to the many-body en ergy 
density was estimated in the low density limit [2(| l2lf . 
It was further discovered that the energy density of uni- 
tary Fermi gases can be directly related to an asymptotic 
property of the reduced two-atom correlation function, 
which illustrates the important role of universal few-body 
physics in the many-body limit [2^ - [25| . The interest- 
ing interplay between two-body physics and many-body 
physics explored in spectroscopy experiments was em- 
phasized 2(1 27 1 . A few separated efforts have also been 
made with an emphasis on the other side of the coin, i.e. 
how few-body structures respond to the presence of a 
quantum background, i.e. either statistical or dynamical 
effects. These studies have led to the findings of anoma- 
lous dimers with tunable masses and new universal be- 
haviours of Efimov states in quantum mixtures 28- 3ljj . 
It remained to be seen whether these predictions can be 
tested in future experiments. Furthermore, it was also 
shown that these anomalous few-body states can quali- 
tatively modify the low energy multiple scattering phase 
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shifts 

The approach outlined in this article is an alternative 
to the more elaborated diagrammatic rcsummation that 
was previously carried out by the authors and their col- 
laborators in Ref. 33j . The two approaches yield almost 
identical results. In section II, we first carry out a sim- 
ple scaling argument, as a caricature of resonating Bose 
gases, illustrating the relevance of fermionization in Bose 
gases near resonance and discuss the limitation of the 
coarse grain procedure. We also remark on a close re- 
lation between the effective two-body interactions and 
Lee-Huang- Yang corrections in Ref. [6]. In section III, 
we discuss the energetics of dimers and trimers in con- 
densates and explore the implications on the many-body 
physics. Especially, we point out that in addition to the 
fermionization phenomenon, an instability sets in at a 
positive critical scattering length as a signature of forma- 
tion of dimers in condensates. We also show this aspect 
is a consequence of the sign-change of the renormalized 
two-body interactions between the condensed atoms; the 
effect of condensates on the two-body interaction con- 
stant is investigated via taking into account the self en- 
ergy of dimers and via imposing an infrared boundary 
condition for the renormalization flow. In section V, we 
summarize the results of diagrammatic resummation. We 
employ a set of self-consistent equations to address the 
two issues listed at the beginning of the article. These 
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self-consistent equations were first introduced by the au- 
thors and their collaborators in Ref. 331. To make it 



accessible to general readers, we leave the most technical 
part in the appendix and outline the framework of the 
approach and result in this section. In section VI, we 
conclude our studies. 



atoms interact in a box of size £ at arbitrary scattering 
length a. £2(£, a; 0) for a contact resonance interaction 
can be worked out, either by assuming two atoms are in 
a harmonic trap of harmonic length £ or in a block of 
size £. The asymptotic behaviours are universal up to 
numerical prefactors. For two atoms in a block of size L, 



A CARICATURE 
Relevance of Fermionization: a scaling argument 

The energetics of Bose gases near resonance can be 
qualitatively understood via a coarse grain procedure 
which is more or less equivalent to the real space renor- 
malization transformation. The simplest implementation 
of this is to first divide a quantum gas into N blocks each 
of which is of the size of £ x £ x £ where £ = 1/ y/2mjl 
is the coherence length and /i is the chemical potential. 
Because the chemical potential is a non-additive thermo- 
dynamic quantity, it is natural to define it as the change 
of energy when adding an additional atom to a particu- 
lar block and the effect of other blocks is to set an ap- 
propriate boundary condition. Therefore, the chemical 
potential can be considered as the interaction energy be- 
tween the added atom and existing atoms in the block. 
If we further assume the interaction energy is dictated by 
a pairwise one, then 



A* = £2(£,a;n)n£ 3 



(1) 



where e2(L = £, a; n) is the characteristic interaction en- 
ergy between two atoms in the block and n£ 3 is the num- 
ber of atoms in the block with n being the number den- 
sity. This is a standard coarse grain procedure which 
relates a microscopic quantity e2(£,a;rt) and a thermo- 
dynamic quantity, the chemical potential fi. The esti- 
mate of €2(£, a; n) itself is a full many-body problem that 
is usually very difficult to carry out. In the dilute limit 
however, one can show that when two atoms interact in 
a box of size £, the probability of being scattered by the 
third particle is negligible because the mean free path I 
is proportional to l/n47ra 2 (a is the scattering length) 
much longer than £. In fact, 



(2) 



which is small in the low density dilute limit. So at least 
in this limit, we can approximate t2(L — £, a; n) as the 
energy of two interacting atoms e^{L = £, a; n — 0) in 
an empty box of the size of the block. If we assume 
this is also qualitatively correct even in the unitary limit, 
then we have a very simple self-consistent equation; the 
only knowledge we need to solve this equation is how two 
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when a -C i; 
when a> L. 



(3) 



C± t 2 are two positive prefactors depending on the details 
of the block and arc of little importance for our qualita- 
tive discussions here. It is important to notice that at 
resonance, e 2 (£ = £, a = oo;0) is finite and scales like 
the kinetic energy of an atom moving in an empty box of 
size £. 

Substituting the results in Eq. ([3]) into Eq. (JlJ, one 
obtains the estimate of chemical potential. In the dilute 
limit, 
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(4) 



the correction to the first Hartree-Fock term is the lead- 
ing finite size correction to the interaction energy and 
belongs to the well-known Lee-Huang- Yang effect. When 
a approaches infinity on the other hand, this simple pro- 
cedure leads to the prediction of fermionization. That 
is 
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(5) 



which scales as the Fermi energy of a quantum gas with 
the same density and mass. Although crude, the coarse 
grain shown here points to a phenomenon that was pre- 
viously seen in a few numerical calculations. Given that 
it is very simple, we consider it quite a success. The rele- 
vance of fermionization to Bose gases near resonance was 
observed in a few theoretical studies 37-4oj|. 

This aspect of Bose gases near resonance is also an 
essential feature of Tonks- Girardeau gases or hardcore 
bosons in one dimension 3J, 1 3 -51] . The one-dimensional 
Bose fluids were later further studied using the Luttinger 
liquid formulation |36j; the renormalization of micro- 
scopic parameters also plays a very important role in the 
construction of that theory although the implementation 
was carried out very differently from what we are going 
to focus on below. 



Running two-body interaction constants 

But how good is the starting point that near resonance 
we can approximate £2(X,a;n) as the two-particle inter- 
action in an empty box completely neglecting the effect 



of many other identical particles? To address this, we es 



timate e 2 (L, a; n), the interaction of two atoms in a bo: 
of size L via employing a more sophisticated approach 
the real-space renormalization transformation (RSRT 
which further takes into account the many-body effec 
on e 2 {L, a;n). This approach indicates the fermioniza 
tion must not be the whole story. 

Consider, instead of e 2 (L,a;0) discussed above, 




g 2 {L,a;n) = e 2 (L,a;n)L 3 (6, 

which is the effective strength of the short range two- 
body interaction. Again we divide the length scales in 
RSRT into two regimes that are separated by £ : the short 
distance regime in which two- and few-body physics dom- 
inates and the long wave length regime where the many- 
body collective effect dominates. £ defines the interface 
where the microscopic few-body parameter g 2 at shorter 
distance needs to match the macroscopic coarse grain 
condition. 

So at scales smaller than £, we can employ the RSRT 
of the two-body running coupling constant g 2 (L,a;0) 
in vacuum to monitor the effective interaction; the fi- 
nite density has very little effect in this regime i.e. 
g 2 (L,a;n) — g 2 (L, a; 0) + 0(L/£). At larger distances, 
because 172 defined here is subject to a thermodynamic 
constraint of the chemical potential at L ~ £, the effect of 
condensate on g 2 or e 2 is to impose a boundary condition 
on the flow of g 2 (L,a;n) via the coarse-grain condition 
in Eq. (JTJ. And in this approach, the collective physics 
at scales larger than £ influences the flow solely through 
a simple boundary condition. 

Practically, since L defines the size of micro-blocks in 
the renormalization procedure, it therefore defines the 
momentum cut-off via A = L . The transformation 
from L to L' is equivalent to rescale the momentum from 
A to A' = L'^ 1 . To obtain the running coupling constant, 
one can use the standard momentum-shell renormaliza- 
tion procedure to track the transformation from the orig- 
inal g 2 to the new one g' 2 when the Hilbert space or the 
momentum cut-off is rescaled from A' to A = A' — SA 
(See Fig. [JJ). 

The reduced two-body Hamiltonian we use for this il- 
lustration is 



H 2 -body — H < + H > + H >< 

= E t'Mfr' + fi E^-k^-k' + h.c.) 

k k,k' 

H > = E e P b l b p + H E^-pV^-p' + h.c.) 
p pp' 

5 >< = IE^M-p + ^ (7) 



Figure 1. (Color online) a) Schematic of the renormalization 
of the low energy on-shell scattering amplitude. An initial 
state (k, — k) on the inner sphere is first scattered into off- 
shell high energy states (p, — p) represented by the shell region 
between two outer momentum spherical surfaces (dashed) be- 
fore finally being scattered backed into (k', — k') on the inner 
sphere. Here the outer dashed spherical surface is defined by 
A' and the inner dashed one by A; A < |p| < A', b) The 
one-loop diagram that leads to the renormalization equation 
below. The internal lines are for states within the shell re- 
gion described in a). Each vertex stands for the two-body 
interaction g-z (A'). 

Here frj^&k) is the creation (annihilation) operator for a 
Bose atom with momentum k, f2 is the volume. The sum 
in H < is over momenta |k|, |k'| smaller than A and the 
sum in H > is over the states within the shell in Fig. [IJ i.e. 
|p|, |p'| G [A, A']. H >< describes the off-shell scattering 
from low momentum states with |k| < A into the high 
energy states with momenta p within the shell and vice 
versa. This interaction can induce an effective scatter- 
ing between low momentum states (k, — k) and (k', — k'), 
|k|, |k'| < A, via virtual states (p, — p) within the shell. 

When rescaling, the states within the shell thus lead 
to an additional contribution to the two-body interac- 
tion in H < . One can obtain the beta- function for the 
renormalization equation diagrammatically. This calcu- 
lation is very similar to the T-matrix calculation except 
one should restrict to the virtual states within the shell 
between A' and A = A' — SA. The diagram in Fig. Q]rep- 
resents such an additional contribution to the effective 
two-body interaction, 

-tto(A')=.V(A') / (^3/ ^G°(e,p)G°(-e,-p) 



Here the momentum integral J is over the states within 
the shell shown in Fig. [TJ i.e. A' > |p| > A. 

After carrying out the energy and momentum inte- 
grals, one can easily find the transformation for the two- 
body interaction constant 
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g 2 (A> - 5A) = g 2 (A') - ^gli^SA + 0(<5A 2 ). (9) 

The transformation of g 2 under the real space rescaling 
can be obtained by converting A to L _1 . For a quantum 
gas with a finite density, we therefore have (g 2 = 92 /L) 



dg 2 (L) 



g 2 (L = R*) 



U . „ 



(10) 



The boundary condition at Z/ = £ is exactly the condi- 
tion in Eq. ([T]), i.e. at scale £ the microscopic running 
coupling constant has to match the thermodynamic con- 
straint suggested by fi, assuming the main contribution 
to fi is from the two-body interaction g 2 {L,a;n). At a 
very short distance i?*, the boundary condition is set 
by Uq, the strength of the bare two-body short range 
attractive interaction with range R*. And for the reso- 
nance phenomena we are interested in, R* is always much 
smaller than a. 

By contrast, in a vacuum the coupling constant g 2 
should flow to the value of Aira/m, the standard form 
of the two-body effective interaction, or 



dhiL- 1 



m ~2,T\ 



.g 2 (L = i?*) = ^,.g 2 (i->co,a;n = 0) = ^. (11) 

For a bare attractive two-body interaction with strength 
Uq(< 0) and range R* , the boundary condition at L = 00 
in Eq. (| 1 1 1) establishes a well-known relation between Uq 
and the scattering lengths a. g 2 = Lg 2 (L) as the solution 
to Eq llll can be expressed in terms of a, 



g 2 {L,a;n = 0) = 



A-kcl 1 



m 1 



2a 



(12) 



Obviously, g 2 appears to be repulsive only in the limit 
of long wave length when L ^> a. At short distances 
R* < L < a, 



92{L) 



2n 2 L 



(13) 



is negative and independent of R* or a. Eq. (fT3"|) in- 
dicates a universal form of the two-body running cou- 
pling constant that induces resonance scatterings. This 
crossover from repulsive to attractive interactions hap- 
pens at L* ~ a. 

One can further show that for a repulsive interaction 
that leads to the same zero energy scattering length a, 




Figure 2. (Color online) The solution to the self-consistent 
boundary condition, Eq. (|14jl . The solution is obtained by 
solving Y(x) = x 2 /(n 1/3 a) - (2/n)x s = 87r, where x 2 = 
2m/i/n 2 / 3 . From the top to bottom are Y(x) for a < a cr , 
a = a cr and a > a cr . At a cr , the equation has only one 
solution and above a cr there are no real solutions. 



g 2 also flows toward the value of Aira/m when L is much 
longer than the range of interaction. For instance for a 
hardcore potential with a — R* where R* is the radius 
of hardcore, one obtains the same expression as Eq. (|13[) 
except that the range of L is 2a/iT < L < oo; and not sur- 
prisingly, g 2 in this case is repulsive for arbitrary length 
scales. 

So only in the long wave length limit, the attractive in- 
teraction with positive scattering lengths yields the same 
physics as the repulsive ones even though at short dis- 
tances they are distinctly different. At the zero energy, 
the effective interaction is Aira/m, repulsive as far as a 
is positive disregarding whether the bare interactions are 
repulsive or attractive. This for long has been a common 
belief in the field of cold atom physics. 

As we will see below, this no longer holds near reso- 
nance when the many-body renormalization effects due 
to condensed atoms are further taken into account. The 
reason for this is the low energy window where we can 
approximate the resonance interaction as a repulsive one 
(which is of order 1/ma 2 ) gets so narrow that the effect 
of condensates on the two-body coupling becomes par- 
ticularly pronounced near resonance. 

The renormalization-equation approach had been pre- 
viously applied to analyze the effective field theories for 
few-body scattering phenomena 4l[ . They were also em- 
ployed to identify the coupling constants and quantum 
phases in the field theory models for the lower branch 
unitary gases with a lot of successes 42 - 44j . It was later 
employed to explore the physics of geometric resonances 
and confinement induced scattering phenomena (45| . Our 
application to Bose gases is perhaps another excellent 
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example to demonstrate that the simple and generic ap- 
proach of renormalization can lead to some surprising 
breakthroughs . 

Eq. (ITU|) is a RSRT equation which satisfies Eq. JT]) 
and yields an estimate of £2(£, a; n) or g 2 - The boundary 
condition leads to a self-consistent equation for \i. When 
expressing in terms of a using the solution to Eq. dill) , 
one finds 



Ana 



1 



fi = n- 



1 — -^J2m[ia 



(14) 



Again in the limit where a is much less than £, the equa- 
tion yields the Hartree-Fock energy plus the correction 
of Lee-Huang- Yang character; 



The RSRT suggests an important feature that is ab- 
sent in the simplest coarse grain approach (Eq. (jTJ) , (J3j) ) . 
It turns out that near resonance, there is a substan- 
tial modification of the underlying two-body physics, i.e. 
dimer energetics and therefore the interaction energy be- 
tween condensed atoms £2; it can no longer be justified 
to approximate £2(£, a ! n ) as the interaction energy in an 
empty box or at zero density. In fact as shown below, the 
uplifted dimers (towards condensates) cause an instabil- 
ity of atomic condensates when approaching resonance 
from the molecule side. The emergency of the imaginary 
part of the chemical potential beyond the critical point 
signifies a hybridization between atoms and molecules 
which is missing in the simplest coarse grain argument. 
In the next section, we further elaborate on this fascinat- 
ing aspect of Bose gases near resonance. 



Ancm. As/2 / — ? , ... r , 
(1 + ^Vm 3 + ...). (15) 



DIMERS AND TRIMERS IN A CONDENSATE: 
THE SPECTRUM FLOW 



Another solution of /i scales as 1/a 2 consistent with the 
binding energy of lower branch molecules and we don't 
consider here. In the unitary limit however, the equa- 
tion not only indicates fermionization but also suggests 
a critical point beyond which there are no real solutions 
to the equation. This is most obvious when a is infinity 
and f/2 becomes negative. This property of Eq. (|14j) is 
illustrated in Fig. 2. One can show that at the critical 
point, 



^3^= 1 1/3 cr = 27r 4/3 
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,2/3 



(16) 



and two real solutions merge into a single one. Beyond 
this point, the equation yields to a complex solution to 
the chemical potential. 

I 



How are the dimers or trimers formed in the presence 
of a condensate or of a quantum gas? In the context 
of quantum mixtures, there have been a few attempts 
to answer this question: how are the few-body struc- 
tures affected by the presence of a Fermi surface [28143 ll | ? 
Surprisingly so far little efforts have been made to under- 
stand the dimers and trimers in the presence of a conden- 
sate, partially because the background of a condensate is 
more dynamical compared to that of a Fermi sea. Since 
this plays a critical role in the interplay between few- 
and many-body physics that interests us, here we make 
an effort to estimate the effect. 

It is possible to solve the two-body and three-body 
S-matrices in the presence of the many-body effect due 
to the self-energy. Assuming the self-energy of quasi- 
particles is E, one finds that for two incoming atoms with 
momentum p, — p scattered into p',— p' and with total 
frequency E, 



G 2 (E;p,p') = Uo + Uo 



d 3 q 



(2tt) 3 E-2r,- 2e c 



■iS- 



;G2(£;q,p') 



(17) 



where 77 = E — /1. A diagrammatic representation is given 
in Fig. 3(a). In the dilute limit, r\ = fi = Airan/m; in 
general, E, fj. and rj are unknown and need to be de- 
termined self-consistently later on. For now, we simply 
assume that 77 is a given parameter. 

G*2(i?;p,p') = G 2 {E) as a result of the short range 
interaction and note that when rj = or in vacuum, 
G\(E = 0) = A-Ka/ra (superscript indicates the case 
of vacuum) and G^E) = 47ra/m(l — iay/mE). One can 



then show that in this approximation, 



G 2 {E) = G° 2 (E - 2r?); 



(18) 



the pole of G 2 (E) is shifted from the pole in vacuum by 
2rj. The pole defines the dimer binding energy and so 



2-q 



(19) 
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where eo and e° D are the binding energy of dimers in 
the presence of a condensate and in vacuum respectively. 
At a given positive scattering length, the dimer spectrum 
flows (in the energy space) towards the zero energy where 
the condensate lives as one increases the 77. 

One can also calculate the amplitude of three-body 
scatterings corresponding to the processes described in 
Fig. 3(b). We consider a general case where three in- 
coming momenta are ki = p/2 — q, k2 = p/2 + q 
and k3 = — p, and outgoing ones are k' x = p'/2 — q', 



k' 2 = p'/2 + q' and kg = — p'. The scattering ampli- 
tude between theses states can then given by A 3 (E; p, p') 
where q and q' doesn't enter explicitly; it represents the 
sum of diagrams identical to Fig. 3(b). 

It is more convenient to work with the reduced ampli- 
tude G 3 {E; p) = A 3 (E; p, 0) where p' is already taken to 
be zero. G 3 {E\ p) itself obeys a simple integral equation 
as can be seen by listing the terms in the summation ex- 
plicitly. The diagrams in Fig. 3(b) yield (see Appendix; 
the mass is set to be one, i.e. m = 1), 



]g 3 (E, P ) =-l K (E~3 mP ,0) + - 

4 Z 7T 



dq 



K{E~Z mP ,q)q 2 



+ 



dqdq' 



l a 2 

4y 



K(E~3 V ;p,q)q 2 



-1 



f 3r)-E + i5+-±Q +3r > 
K{E-i mqi q')q 12 



1 



E 



^+3r l -E-lJy^^-l^ + ^-E 



(20) 
(21) 



where K(E — 3i];p, q) is the kernel defined as 



K{E - 377; p, q) = — In — , (22) 

pq p z + q z — pq + or) — E 

The sum of the above infinite series leads to the following 
integral equation of G3 as 



G 3 (E,p) = -2K(E - 3r,;p, 0) + - / dq K{E *™,<l)<? G 3 (E, q). 



§g 2 + 3ry - £ 



(23) 



r 



When 77 — as in vacuum, this equation is identical to an 
integral equation previously obtained in an atom-dimer 
model to study the renormalized three-body forces [47j |. 
Comparing to G^{E,p) in vacuum when 77 = 0, again 
one finds that the energy of a trimer in a condensate ct 



is related to e^, its vacuum value via 



3ry. 



(24) 



What Eq. (|24|) shows is a simple fact of a condensate. 
If all the finite momentum atoms have a mean-field en- 
ergy shift S — fj, with respect to condensed atoms, the 
energy of few-body bound states (with finite k compo- 
nents) experiences the corresponding energy shifts. As 
a consequence, when = 0, we should expect that the 
three-body forces in a condensate should be divergent. 
This was observed numerically in Ref. [33j |; the three- 
body potential is divergent when 3r/ = — e n where e n are 
the Efimov eigen values with n = 1,2,3,.... 

What is the consequences of the spectrum flow or the 



energy shift due to the condensate? The main conse- 
quence is that in a condensate, a dimer crosses the zero 
energy, or the energy of condensed atoms at a positive 
critical scattering length ao or cd =0 when 



277(00) = 



(25) 



D 



where r\ itself is a function of an- By contrast, in vac- 
uum, a dimer crosses the zero energy or the scattering 
threshold at resonance or a = 00. If we simply apply the 
Hartree-Fock approximation S = 2/x = Siran/m, we find 
77 = jj, and 



(1/87T) 



1/3 



(26) 



Beyond this point, one has to take into account the hy- 
bridization between atoms and molecules. The dimer 
formation in condensates was previously studied in a 
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a) 



b) 



Figure 3. (Color online) Diagrams for the calculations of 
dimer and trimer energy in a condensate, a) is for two-atom 
channels. Each solid internal line here is for a Green's func- 
tion of non-condensed atoms with the self-energy effect taken 
into account, G _1 (e, k) = e — ek — S + /i + iS. b) is the loop 
diagrams for the three-atom channel. 



random-phase approximation; those results a re q ualita- 
tively consistent with the picture painted here |48| . Pair- 
ing instability and formation of molecules in the upper 
branch Fermi gas was emphasized in Ref. [49j |. 

Since it is necessary to have molecules below conden- 
sates so that condensed atoms have effective repulsive 
interactions, the penetration of dimers into the conden- 
sate implies a change of the sign of interactions, from 
repulsive to attractive ones that can lead to a potential 
instability. Below we further amplify this aspect. 



SIGN CHANGE OF g 2 : 
A CONSEQUENCE OF SPECTRUM FLOW 

In a condensate, the low energy two-body interaction 
constant is renormalized not only by the virtual scatter- 
ing states as in vacuum but also by the interactions with 
the condensed atoms. The later effect is many-body in 
nature. Below we focus on the particular many-body ef- 
fect related to the self-energy of virtual states and include 
this in the renormalization procedure. The self-energy of 
non-condensed atoms is due to scatterings by condensates 
and depends on the interaction strength and density of 
atoms. 

We now apply a self-consistent renormalization group 
equation (RGE) to investigate this issue. To study the 
coupling constant, we start with an assumption that the 
self-energy and the chemical potential of non-condensed 
particles are already given as £ and /i. The simplest 
Hartree-Fock Green's function for virtual atoms is of the 
form 



G(e,p) = 



Pl 

2m 



iS 



(27) 



We then calculate g^i^,^), the zero energy effective in- 
teraction between condensed particles for a given r\ — 



£ — [i using a very similar procedure as that in sec- 
tion H except G° in Eq.® now should be replaced 
with G defined here. The corresponding renormalization 
group equation (RGE) for the running couple constant 
92 = 92 (A) Ah can be found to be, 



dg 2 (A h ) 




(28) 



where is the dynamical length relevant to the renor- 
malization transformation and depends on the many- 
body parameter rj. When 77 = as in vacuum, g% flows 
to the desired value of Aita/m. With a finite 77, we find 
that 52 runs to the following value 



<7 9 = lim 



92{A h ) _ 4wa 



1 



1 — \J2mr\a 



(29) 



as A/j becomes zero and all the k ^ virtual states are 
included in the renormalization transformation. The re- 
sultant g% is the effective interaction between condensed 
atoms after all non-condensed or virtual states are inte- 
grated out. Eq. (|2"9")l was proposed in Ref. 33j as an 
effective interaction for condensed atoms. This is also 
fully consistent with the RSRT result presented in sec- 
tion II. 

At first sight, the structure of Eq. (|2T)1) appears to be 
very similar to the zero-density expression for the two- 
body running coupling constant 52 (L) in Eq. ([12"]) . How- 
ever, the physical implication is entirely surprising. First 
of all, </2i the effective interaction between condensed 
atoms now depends on k v = \J2mr\\ so it is now a func- 
tion of £ — fi, or the density of the gases reflecting a 
many-body effect. In the dilute limit, 



g ° = —(l + V8Trna 3 + ...); 
m 



(30) 



the first term stands for the Hartree-Fock energy and the 
second one yields the Lee-Huang- Yang type correction to 
the energy density of Bose gases. 

And most importantly, unlike in vacuum where the 
zero energy effective interaction constant A^ajm is al- 
ways positive as far as a, the scattering length remains 
positive, in condensates g% (A = 0) is positive only in 
the dilute limit when ^/mqa ~ y/ na 3 <C 1. When ap- 
proaching the resonance, for a given S — ji, the effective 
interaction between condensed atoms becomes negative 
before a becomes infinity as indicated in Fig. 4. In other 
words, the presence of a condensate completely alters the 
flow of the coupling constant at the low energy limit; it 
changes the sign of the effective interaction constant near 
resonance. 
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The property of Bose gases near resonance is dictated 
by this change of the sign of interactions. In fact as a 
precursor of this, pure atomic condensates lose metasta- 
bility as seen in the following more elaborated discus- 
sion and diagrammatic resummation in section IV. Mi- 
croscopically, the change of sign of g% is correlated with 
and driven by the molecules' entering the condensate. In 
the approximation employed here, the sign change occurs 
exactly when the molecules penetrate into the conden- 
sates at scattering length ajj (see Eq. (|25jl). 

To further determine r\ or /i and E and understand 
the effect of the sign change of gi on the condensates, we 
should specify the RGEs with a boundary condition. The 
following steps have to be carried out. Once gi is found 
as a function of E and /1, one can apply it to calculate 
E(riQ,ii), the energy density of the system with uq con- 
densed atoms and non-condensed particles at chemical 
potential /1. Following the general thermodynamic rela- 
tions, the chemical potential for the condensed particles 
fi c should be 



_ dE(no,fi) _ 1 o 2 



(31) 



And for the ground state we further require that the con- 
densed atoms are in equilibrium with the non-condensed 
reservoir at chemical potential /i. 

/i = u c (32) 

which was first suggested by Pines and Hugenholtz jiij . 
One can verify that Eq. (|29|) .(|31 |) and (|32"T) are identical to 
the corresponding self-consistent diagrammatic equations 



employed in Ref. [33j. More explicitly, one finds that for 



9° n o 



He 



(33) 



dglcm 

drj dno 

One can view Eq. (|3"T|) and (j3"2"| as a boundary condi- 
tion for 32 (Ah) in the RGE in Eq. when Ah = and 
if <9E/<9no is given. To finally solve the equation self- 
consistently, one needs to supply Eq. (|40[) in the next 
section to further determine that d~E/dno = tj/hq and 
the set of equations produced in this way are identical to 
the set in Ref. [33| derived diagrammatically. 

Here to illustrate the main features, we make a few fur- 
ther simplifications without losing the generality. One is 
that we neglect the no-dependence in E so that fi c — 
g^riQ- Second is that we further approximate Uq as n be- 
cause they are of the same order in the regime of our in- 
terest. We then have a single parameter renormalization 
equation Eq. ([28]) with the following boundary condition 



r 92(A h ) fi 
<k = hm — = — . 

Ah^o Ah n 



(34) 



3 

-5 
20 

a? 

-20 



(b) 




1/an 



1/3 



Figure 4. (Color online) a) <72 (in units of 4na/m) as a func- 
tion of k v = ^J2mr\ at a given scattering length a; g — £ — fi 
is a function of density and is equal to iTrna/m in the di- 
lute limit. Note that when ak v = V 8Tvna s <C 1 or in the 
dilute limit, Q2 approaches its vacuum value of 4na/m but 
deviates from it substantially once k v a is of order of unity. 
At resonance when 1/a — 0, §2 is negative for any arbitrary 
r\ implying attractive interactions between condensed atoms, 
b) Illustration of the dimer energy (in units of n 2 ^ 3 /2m) in 
the presence of a condensate (the upper curve). The dashed 
line indicates dimers are no longer well defined because of the 
coupling to the continuum. As a reference we also show the 
dimer energy in vacuum (the lower curve). Note in vacuum, 
the dimers reach zero energy right at resonance. 



Last, although 77 = E — /i in general should be /3/j, with 
(3 being an unknown but smooth function of a, no and 
/i, in the dilute limit (3 = 1. Eq. (|4T))) below implies 
that j3 varies between 1 in the dilute limit and 2/3 in 
the fremionized limit that interests us. So we can neglect 
its variation by simply setting f3 — 1 for this part of 
discussion. Eq. (|2"5|) and the boundary condition for the 
RGE in Eq. (f3~4"]) now lead the following single parameter 
self-consistent equation for /x, 



n 



Ana 



ni(l — y/2mfla) 



(35) 



which, apart from a numerical prefactor in the denomi- 
nator, is identical to Eq. (I14[) which was obtained empir- 
ically. The numerical solution of this is presented in Fig. 
5. 

Two essential features are shown in Fig. 5. First, 
the chemical potential reaches a maximum at 
a precursor of the sign-change of two-body interaction 
32 near resonance. The value of maximum is around 
89%e^ , very close to the values obtained in a constrained 
variational approa ch [38l| and in a diagrammatic resum- 
mation approach [33j (see Table I for details). Here 
£f = (6tt 2 ) ' 
the number density n 



2N > 2 / 3 n 2 / 3 /2m is the Fermi energy for a gas with 
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1 



3. 




Figure 5. (Color online) The numerical solution to the self- 
consistent equation Eq. (|35[) . The chemical potential (the real 
part) reaches the maximum (blue circle) when n 1/ ^ 3 a = 0.18; 
beyond this point the chemical potential develops an imag- 
inary part (dotted line). The dashed line is the chemical 
potential in the Lee-Huang- Yang theory. The smooth con- 
tributions from the three-body potential gz (not shown here) 
were studied in the previous diagrammatic calculations [33J 
and turn out to be around a few percent of the effect shown 
here. 

And above the critical scattering length, the chemical 
potential develops an imaginary part. 





Fermionized* 


LHY 


Efimov Max. in fi, 






effect physics 


instability 


Cowell et al, 


Yes 


No 


No 


No 


2002 [37] 


(2.92ejO f 








Song et al., 


Yes 


No 


No 


No 


2009 [38] 


(0.80ef) 








Lee et al., 


Yes 


No 1 " 1 " 


No 


No 


2010[39] 


(0.66e F ) 








Diederix et al., 


Yes 


No n 


No 


No 


20 11 [40] 


(0.83£f) 








Borzov et al., 


Yes 


Yes 


Yes 


Yes m 


2012[33J 


(0.93£f) 









Table I. Comparison of Different Theory Approaches 



* The lower bound of /j, was measured to be around 
0.44e F in the ENS experiment [11]. 

| The value in the bracket indicates the estimated chemical 
potential; same below. The estimated chemical potential 
2.92ef exceeds the result for a completely fermionized gas. 
ff In the field theory approaches there, the signs of the 
correction of order of V na 3 are opposite to the LHY effect. 
However, the LHY effect was reproduced in the numerical 
program in Ref . [40l] . 

f f f This is seen both in the diagrammatic resummation and 
the RG approach outlined here. Note that 0.93eF is for a 
range of three-body parameters relevant to cold atoms. 



Imfi 



8 ^ ^ a 
(37r) 2 / 3 a cr 



1) 



1/2 



(36) 



when the scattering length a is increased slightly beyond 
the critical point a cr = 0.18n -1 / 3 . 

A different renormalization group approach based in 
an atom-molecule model was also applied in a previous 
study to understand Bose gases near resonance [391 ] . Our 
results differ from theirs in two aspects. First, in our ap- 
proach, an onset instability sets in near resonance even 
when the scattering length is positive, a feature that is 
absent in that previous study. Second, when extrapo- 
lated to the limit of small na 3 , the results in Ref. [39j 
imply a correction of order of Vna 3 to the usual Hartree- 
Fock chemical potential but with a negative sign, oppo- 
site to the sign of LHY corrections and/or our results. 
The results of the self-consistent approach in Ref. 
are similar to the ones in Ref. [39[ but differ from ours. 
In table I, we make further comparisons by listing the 
main features in different approaches. 



DIAGRAMMATIC RESUMMATION: A 
SELF-CONSISTENT APPROACH 

From a phenomenological point of view, it is quite ap- 
pealing to generalize the self-consistent coarse grain re- 



lation Eq. ([T]) by further taking into account the three- 
body effective interaction 53; 



l« = ng 2 (£, a; n) + —g^, a; n) 



(37) 



where g 2 ^{L, a; n) are the renormalized two- and three- 
body interaction constants respectively at length scale L 
and £ _1 = y/2mjl. If one can calculate these renormal- 
ized quantities, then one is able to obtain \x which in- 
cludes the effect of 53. We have proceeded further from 
here using the renormalization group equations similar 
to what was discussed in section II, III; they yield qual- 
itatively the same results as the diagarmmatic approach 
below. However, when benchmarking against the dilute 
gas theory, the diagrammatics turn out to be numerically 
superior; the diagrammatic resummation used in our pre- 
vious paper reproduces 99.96% of Lee-Huang- Yang cor- 
rections in the dilute limit. For this reason, we are going 
to outline the framework of the diagrammatic calcula- 
tions and briefly comment on the results. 

The Hamiltonian for a condensate with a short range 
interaction is 
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H = ^(ek - n)blbk + 2U n ^ &[&k 

k k 



2^/VL 



k',q 

t fct 



h.c. 



+ 2^ Z 6 lc+a &T -k+! 6 k'+l & -k' + a+^C- (38) 

k,k',q 

Here ek = k 2 /2m and the sum is over non-zero mo- 
mentum states. Uq is the strength of the contact in- 
teraction related to the scattering length a via Uq 1 = 
m(47ra) _1 — 5Zk(2ek) _1 , fi is the volume, no is the 
number density of the condensed atoms and /j, is the 
chemical potential, both of which are functions of a and 
to be determined self-consistently. 

In the diagrammatic approach, we first define the 
chemical potential of non-condensed particles as /_t and 
the number density of condensed atoms is Uq. The en- 
ergy density can be calculated as -E(n ,u) (see below). 
Then one should have the following set of self-consistent 
equations for a gas with total number density n, 



dE{n ,n) 
dn 



n = no — 



A* = fJ>c 



(39) 



where fi c is the chemical potential for the condensed 
atoms and has to be equal to the chemical potential (i 

I 



in equilibrium. We further introduce the self-energy £ of 
non-condensed atoms or virtual particles to facilitate the 
calculation of E(rio,u) that now explicitly depends on 
E(no, u). Thus, Eq. (|39j) has to be further supplemented 
by 



E(no,/i) = Hc{no,n) 



du c 
d In Uq 



(40) 



which can be proven in the same fashion as the Pines- 
Hugenholtz theorem [Zfjj]. Eq. ((39]) and ((40]) have been 
applied to obtain the chemical potential in 3D Bose gases 
near resonance j33|. 

Calculations of E(uq, /i) for a given £(= r\ + u) can be 
carried out diagrammatically. If we restrict ourselves to 
the virtual processes involving only two or three excited 
atoms and truncate the Hilbert space accordingly, then 
diagrammatically we only need to collect the diagrams 
which contribute to the effective two- and three-body in- 
teraction constants 52,3- As far as the chemical potential 
is concerned, this truncation turns out to be highly pre- 
cise in the dilute limit. The result is listed below. The 
mass to is set to be one. 



E(no,fi) = ^77^2(277) + -^n 3 g 3 (3Tj) 

33 (377) =6 ff2 2 (3r,)i?e- / dq K ^^ G ' 3 (-3r,,q) 

(41) 

where G 3 (—3r],p) is a solution of the following integral 
equation 



G 3 (-3 V ,p) = - 

7T 



dq 



K(-3r];p,q)q 2 



377 • 



-1 



1 q 2 + 2r\ 



G 3 (-&i,q)] 



(42) 



And — 1/2K{— 3rj;p, q) is again the one-particle Green's 
function projected to the S'-wave channel; it is defined as 



X(-3r ? ; P)g ) = lln P + g (43) 
pq p + q + orj — pq 

In the appendix, we outline the major steps in the loop 
summation which leads to E(no, fJt). 

The numerical solution to these self-consistent equa- 
tions was shown in Ref. [33| and they are qualitatively 
the same as the solution to the self-consistent RGE for 
t/2 and we are not going to repeat here [501 ] . Here we 



I 

want to make a few further comments on the resumma- 
tion technique. 

First, <?2 defined this way is an effective two-body inter- 
action rcnormalized by condensates and includes a subset 
of A^-body interactions defined in the vacuum. At one- 
loop level, it yields the most dominating contribution; the 
residue effects are from the irreducible N = 4,6, ...-body 
interactions which contains less than one thousandth of 
the total contribution. 

Second, the three-body contribution in our self- 
consistent approach appears to be around a few percent 
and numerically insignificant. Since when compared to 
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<72i the contribution from 33 in the dilute limit as well 
as near resonance is small as shown previously 33[, it 



is reasonable to conjecture that further inclusion 54.5,. 
wouldn't change our result presented here in a substan- 
tial way. The truncation of the energy density expression 
at <?3 should be accurate enough for all the practical pur- 
poses of studying Bose gases near resonance. We hope 
these statements can be tested in precision measurements 
of chemical potentials as well as in future quantum Monte 
Carlo simulations. 

Third, the energy density expression in Eq. (I4ip be- 
comes exact in the limit where only the processes involv- 
ing two or three virtual atoms are allowed. Effectively, 
this is equivalent to truncating the Hilbert space and in- 
cluding the correlations up to the trimer channel. 



CONCLUSIONS 

The RGE approach is instrumental to our understand- 
ing of the emergent phenomena in quantum few- and 
many-body systems. The application to Bose gases near 
resonance perhaps is another example of what a simple 
RGE transformation can lead to. We have applied this 
approach to understand the nature of Bose gases near 
resonance and found that energetically, the Bose gases 
close to unitarity are nearly fermionized before an on- 
set instability sets in, i. e. the chemical potentials of the 
Bose gases approach that of the Fermi energy of a Fermi 
gas with equal mass and density. Beyond the instabil- 
ity point, the chemical potential has an imaginary part 
indicating strong hybridization with molecules. 

The model we have employed to study the Bose gases 
near resonance is a short range attractive potential which 
has a range much shorter than the interatomic distance 
of the gases or effectively a contact potential. This is 
a very good approximation of real physical interactions 
between cold atoms. If the potential is a short range 
but repulsive one, then Bose gases are always in the di- 
lute limit because the scattering lengths are bounded by 
the range of interactions, disregarding the strength of 
potential. For bosons interacting with a repulsive po- 
tential but with a range comparable to the interparticle 
distance, we should anticipate the physics in this limit to 
be very similar to what happens in liquid 4 He (HI. [51- 55 1. 
The excitation spectrum should develop roton minima 
that imply strong short range crystal correlations. When 
the range of interactions is further increased, eventually 
there should be a quantum transition to a crystal where 
all bosons are depleted from the condensate. The physics 
of repulsive bosons and liquid He belong to a different 
universality class which fundamentally differs from what 
we described in this article, i.e. the properties of nearly 
fermionized Bose gases near resonances with a contact 
interaction. 
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APPENDIX: TWO- AND THREE-BODY SCATTERING AMPLITUDES IN A CONDENSATE AND 32,3 

In the following, we show the explicit calculation of the three-atom scattering amplitude G3, and E(jiq,/j,) for a 
given E(= rj + /1) by adding the diagrams with minimum number of virtual particles involved. The model Hamiltonian 
could be written as: 

H = ^(e k - /i)6] c 6 k + 2U n ^ b\p^ + hj n ^ &t&l k + hj n &k&-k 

k k k k 

+ ^V^J2 6 q & k<+a6_ k , +§ + h.c. + kk+a^k+akk'+ffe-k'+a + h.c, (44) 

where the sum is over non-zero momentum states. Uq is the strength of the contact interaction which is related to 
scattering length a as: 

1 _ m j_V^ _j_ (Ar\ 

(To ~ 4^ ~ flV2ffc' 

k 

where f2 is the volume. Taking into account only two and three body interactions, the energy density could be written 
as: 
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E(n , M) = lnlg 2 (2r]) + ^n^ 3 (3r?), (46) 

where g 2 and 33 are irreducible 2 and 3-body potentials respectively. g 2 could be found by writing the Bethe-Salpeter 
equation as: 

g 2 (2r 1 )- 1 =U - 1 - i J ^^ G (o;,k)G(-a;,-k), (47) 

where G(cj, k) _1 = lo — — n + iS + is the interacting Green's function. So, the two-body potential could be obtained 
as: 

*™ - j^mtr, < 48) 

Similarly, g 3 could be estimated by summing up all N-loop diagrams with 3 incoming and outgoing lines which are 
depicted in Fig. 3(b). We consider a general case where three incoming momenta are ki = p/2 — q, k 2 = p/2 + q 
and k 3 = p, and outgoing ones are k^ = p'/2 — q', k' 2 = p'/2 + q' and k 3 = — p'. The scattering amplitude between 
theses states is then given by A 3 (E; p, p'). At the tree level, the effective 3-particle interaction is: 

r(°) = 1° . (49) 

E - uj m - uj out - 6 p+p < - 77 + id+ ' 

where Wi„ and uj out are frequencies of lines with momentum k 3 and k 3 respectively. Furthermore, fo is the product 
of the perturbation factor f p , vertex factor f v and symmetry factor f s which will be explained later. 

To keep the notations simple, we set to = 1 from now on. Wc consider on-shell limit and substitute Wj„ and u out 
by p 2 /2 — r\ and p' 2 /2 — r\ respectively. To project into the S'-wave channel, we take the average over all directions, 

p(0) -fa, ( p 2 +p' 2 +pp' + 3r]-E \ -fo v(1 , o n /. n \ 

r = Tr~i ln 2 1 ,2 rT~q ? ] = -^-K(E -3n;p,p), (50) 

2pp \p z +p' A -pp' + iri-E) 2 



where wc have defined the kernel K as: 



K(E - W) - ~,ln P 2 2 +Pl+PP' + ^~Z ] (51) 
VP \p 2 +p' 2 -pp' + 3r/- E ' 



Perturbation factor, f pi comes from the expansion of the exponential of interaction term in Hamiltonian in pertur- 
bation method. The diagrams with I vertices could be written as: 

\{V A + V B )\ (52) 

where Va and Vb stand for interaction terms corresponding to different types of vertices namely A and B. For 
example, for a diagram with 2 vertices of type A and one vertex of type B, f p is the factor in front of V^Vb term in 
numerator of the above equation divided by l\. In the vacuum case, where all the vertices are the same and all the 
lines can have non-zero momenta, f p is simply equal to (1//!). 

Vertex factor, f v , is defined as the product of the factors in front of g 2 for different vertices shown in Hamiltonian. 
For vacuum case, all the vertices have 1/2 factor and /„ is equal to (1/2)'. 

The last factor is symmetry factor, f s which shows the number of identical diagrams generated for a given number 
of vertices. For vacuum case, f s — l\ x A 1 where l\ shows the number of permutations of vertices and 4 Z is the number 
of different ways of connecting vertices (2 for incoming lines and 2 for outgoing lines). 

So, in general the prefactor appearing in r'™' where n is the number of loops (n — I — 2) would be 

./•„ = 2("+ 2 >; (53) 

and f = 4. 

could be written in terms of kernel defined above as: 



/d 3 k 1 k 2 1 

— {--K{E - 3n;p, k))g 2 {E - 3n - —;k){--K(E - 3t?; k,p% (54) 
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where the integral over internal frequency has taken, g 2 (w; Q) has the following form in 3D: 

g 2 (uj;Q) = - j (55) 

where above equation reduces to Eq. (|48p in the limit of zero energy and momentum. The effective three-body 
interaction then could be obtained by summing over r^™^s: 

r eff = -2K(E-3r);p,p') + - f dkk 2 K ( E - 3r r>P> k ) x K{E - 3r);k,p') 

dkk* [dk>k» K & x xK(E 3rj;k',p>), (56) 

4 



The sum of above infinite scries leads to the following integral equation for scattering amplitude A3 : 

A 3 (E;p,p') = -2K(E - 3 m p,p') - - [ dkk 2 K{E Z^M ^=A 3 (E- k,p') (57) 

One then obtains the reduced scattering amplitude G^{E; p) = A 3 (E; p, 0) in Eq. (|23)) . 

For the calculation of g$ for condensates, one has to exclude the tree level diagram that no longer exists because of 
momentum conservation. The sum of the rest infinite series leads to the following integral equation for the scattering 
amplitude A^{E] P ,p'), 

A s (E;p,p')=* f dk ^- 3 ^ fc ) fc2 (lK(E-3r,;k,p')-A s (E;k,p')) (58) 



4 



E + 3tj 



The above scattering amplitude could then be applied to calculate the scatterings between condensed atoms when 
setting E, p and p' to be zero in the above equation but with two further modifications. The first change is the 
numerical factor in front of the effective interactions. This factor in condensate case is 1/4 of the factor in vacuum 
case, because there is a 2 x 2 factor for changing the external legs of the external vertices for non-zero incoming 
momenta. So we will get the same integral equation, but the first term in the bracket of the integrand of Eq. (|58|) 
would be substituted with 1/2K. 

The second change would be in the shift of the energy. If we set the momentum of external legs to zero from the 
beginning, In the on-shell limit there is no shift of the energy for uji n and Lu out in Eq. (|49l) . So, the energy of the first 
and the last kernel in all the terms of the Eq. ((55)) other than the tree level term would be E — 2i] in condensate 
case. Subtracting the one-loop contribution which has already been counted in the renormalized g 2 and taking into 
account the above two modifications, we obtain g 3 via the relation shown in Eq. (39), (40) [Hof . 

I 



